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ABSTRACT 

As part of the OzDES spectroscopic survey we are carrying out a large scale reverber- 
ation mapping study of ^500 quasars over five years in the 30 deg 2 area of the Dark 
Energy Survey (DES) supernova fields. These quasars have redshifts ranging up to 4 
and have apparent AB magnitudes between 16.8 < r < 22.5 mag. The aim of the sur¬ 
vey is to measure time lags between fluctuations in the quasar continuum and broad 
emission line fluxes of individual objects in order to measure black hole masses for a 
broad range of AGN and constrain the radius-luminosity (R — L ) relationship. Here we 
investigate the expected efficiency of the OzDES reverberation mapping campaign and 
its possible extensions. We expect to recover lags for ~35-45% of the quasars. AGN 
with shorter lags and greater variability are more likely to yield a lag, and objects 
with lags <6 months or ~1 year are expected be recovered the most accurately. The 
baseline OzDES reverberation mapping campaign is predicted to produce an unbiased 
measurement of the R — L relationship parameters for H/3, MgnA2798, and CIVA1549. 
However, extending the baseline survey by either increasing the spectroscopic cadence, 
extending the survey season, or improving the emission line flux measurement accu¬ 
racy will significantly improve the R — L parameter constraints for all broad emission 
lines. 

Key words: Galaxy: quasars: general - Galaxy: quasars: supermassive black holes - 
Galaxy: quasars: active - Cosmology: dark energy 


1 INTRODUCTION 

There is good evidence that supermassive black holes 
(SMBHs) are present at the centre of all massive galax¬ 
ies (e.g . iKormendy fe Richstonel 1 19951 : iRichstonel Il998l : 
IFerrarese fe Fordl 120051 ). and that there are tight, empir¬ 
ical relationships between the mass of the SMBH and 
properties of the host galaxy , such as stellar velocity 


properties ot tne h ost galaxy, such as stell ar velocity 
dispersio n (IFerrarese fc Merrittll2000l : iGebhardt et al.ll2000l; 
iMcConnell fe Ma 20131) . light concentration ( Graham et alJ 


20011 ). and bulge luminosity an d stell ar mass (jRichstonel 
_|Kormend Y_£^ l200ll : iMarconi fc Huntl 


IMcConnell fc Ma 20131 ). These relationships sug¬ 


gest an interplay between black hole growth and galaxy 
evolution; however, the true nature of this relationship 
is still unknown and is a major ar ea of resea r ch in 
understand ing gal axy e vo lution (e.g. Ikin 3 :2011:1:. 2001: 
Ferrarese fc Ford 2005ff_ Murray, Quataert fc Thompsonl 


2005|_ _ Di Matteo. Springel fc Hernauistl 20051 : 

Di M atteo et al.l 120081 : iPark et al.l 120151 ). Tb better com¬ 
prehend the origin and evolution of the SMBH-galaxy 
relationship and the growth of SMBHs over cosmic time, it 
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is necessary to obtain accurate and precise measurements 
of black hole masses. Direct measurements of black hole 
masses through stellar or gas dynamics require high spatial 
resolution and are therefore limited to the local universe. 

Active galactic nuclei (AGN) provide an alternative 
method of black hole mass measurement. Continuum emis¬ 
sion from the accretion disk is absorbed by gas deep within 
the gravitational potential of the black hole. The broad 
line region (BLR) gas reprocesses this radiation and emits 
Doppler broadened emission lines. The emission line lumi¬ 
nosity varies in response to changes in the continuum emis¬ 
sion in a roughly linear fashion with an associated time lag, 
t, which is the mean light travel time from the accretion 
disk to the BLR, at the responsivity weighted mean dis¬ 
tance, R = ct. The measurement of this time lag, through 
detailed comparison of the emission line and continuum flux 
variations, is referred to as reverberatio n mapping (RM, 
iBlandford fc McKeelll982l ; |Petersonlll993l L 

If the BLR gas is in virial equilibrium and its motion 
is dominated by the gravity of the SMBH, the mass of the 
black hole is 


A/bh 


fcrAV 2 

G 


(1) 


where G is the gravitational constant, r is the measured re¬ 
verberation time lag, AV is the line of sight velocity disper¬ 
sion of the BLR gas estimated from the emission line width 
of the RMS variance spectrum, and / is a dimensionless 
virial factor that converts the measured line-of-sight virial 
product into the true black hole mass. The virial factor de¬ 
pends on the geometry, kinematics, and orientation of the 
BLR, and although it is of order unity, it is expected to differ 
between quasars. 


Reverberation map ping has yielded lags for approx 


imately 50 AGN (e.g., Petgrsgi^etjilJ 


2009al ; iDennev et al.ll2Q10l ; iGrier et al.l 


20041: 


20121 : 


Bentz et al 


Bentz et al 


The lags exhibit a tight power l aw relation¬ 
ship with the continuum lum inosity, XL\ dKaspi et al.l 
l2000l : iBentz et alj l2009al . l2013h . as predicted from sim- 
ple single-photon phot oionization physics dDavid son 1972; 
iKrolik fc McKeel ll978l L This strong correlation is the 
basis for single-epoch black hole mass esti mates based 
on a single epoch of spectroscopy (e.g. lLaoij Il998 


WajideLPetersonfc^falkanl 19991 : iMcLure fc Jarvia 


2002 


Vestergaard fc Petersonl 2006f l. and enables black hole 


masses to be estimated for a far larger sample of AGN than 


Vestergaard 20041: 

Vestergaard & Osmei) 120091: IShen et al.l 

201 ll: Kellv & Shen 

2013)1. 


The single-epoch mass estimation t echnique has 
been used widely i n cosm ol ogy (e.g., ^jUgttejtjilJ 


2010 :: jMortlocl^etalJ [201ll : Vestergaard ct al 


200S 

201C 


Vestergaard fc Osmei) 2 0091: ISc hulze fc Wisotzkil _ 

Trump et al J l201ll: IShen fc Liul 120121 : iKellv fc Shenl 1201.4) . 

As a consequence, it is important to determine the R — L 
relationship accurately and precisely, and any dependencies, 
because both random and systematic uncertainties in the 
R — L relationship are necessarily transferred into the 
single-epoch masses and all subsequent studies. The origin 
of the scatter around the R — L relationship ha s been 
investigated by sev e ral authors dBentz eh_alJ l2009al . 1201.4 
IWatson et al.l l201ll : iKilerci Eser et alT 2015l i , and recent 


evidence suggests an additional depen dence on the R — L 
relationship with the Eddington ratio jDu et al.ll2015h . 

Despite ‘dark energy’ appearing to be t he dominant 
erg y c omponent o f the universe le.g. IConlev et al ' 


2011 

2012 


Blake et al. 


Anderson et al. 201 


Padrnnnabhanet_al 

iHinshaw et al.l 120131 : 
Planck Gollaboration et al.l 12014 1. its nature is unknown 


and gaining additional understanding of its properties 
remains a high priority. It may be possible to reverse the 
R — L relationship and use the time lag to infer the intrinsic 
luminosity of an AG N, and therefore its luminosity distance 
JWatson et al.ll201lh . The resulting distance measurements 
can be used to independently probe the acceleration of 
the universe and dark energy. The high luminosity and 
prevalence of AGN would make them valuable probes 
of the expansion history of the universe over a greater 
redshift range than other methods, and thus could be 
valuable for investigating the time ev olution of dark energy 
dCzernv et al.|[2013 : iKing et al.ll2014l l. 

Previous reverberation mapping campaigns have gen¬ 
erally only observed a small number of AGN in a single 
campaign (~10), using small telescopes, and the majority of 


(e.g.. Clavel et al.l 

199ll; Robinson! 1994 1Wanders et all 

19971: ICollier et al 

19981: IPeterson et al.l 19981. 

19991: 

Kasoi et al. 

200(4 

D eterson et al. 

2002; 

Kaspi et al. 

20071; 

Bentz et al. 

2009b!; IDennev et al. 

2010; 

Barth et al. 

201l|; 

Rafter et al. 201 ll: 

Grier et al. 20121; Rafter et al. 

20131: 

Du et al. 2014ll. As a consequence, thev have focused on 


the brightest and most variable objects, leading to a bias 
towards local, low-luminosity AGN. While there have also 


been long, multi-year c a mpaigns (up to 8 ye ars Peterson 
iKaspi et all 120001 : IPeterson et al.l 120021 : iKaspi et al 


1999 


200711 . they only monitored a small number of quasars so 


the bias towards the brightest and most variable objects re¬ 
mains. Both higher redshift and higher luminosity quasars 
have longer lags due to time dilation and the R — L relation, 
and higher luminosity quasars also have lower variabilit y 
amplitudes (IVanden Berk et al.ll2004l : lMacLeod et ahlboidl 'l. 
To account for these effects we require larger telescopes and 
longer observation campaigns. This Chapter investigates a 
large scale reverberation mapping campaign being run as 
part of the ongoing Dark Energy Survey (DES), in con¬ 
junction with the OzDES spectroscopic survey. This 5 year 
campaign covers a large range in magnitude and redshift, 
allowing reverberation mapping studies of a much broader 
AGN sample over a large redshift range. 

DES is an optical survey aimed at understanding the 
expansion of the universe using four complementary meth¬ 
ods: Type la supernovae (SNe), baryon acoustic oscillations, 
weak lensing, and galaxy cluster counts. DES officially be¬ 
gan in the 2nd half of 2013, and it plans to image 5000 
square degrees with 5 filters (g,r,i, z,Y) over 5 years. The 
supernova component of the survey will consist of repeated 
observations of 30 square degrees of sky in the g, r, i, z fil¬ 
ters, divided into two deep and eight shallow SNe fields, to 
detect and monitor supernova and other transients. 

OzDES is the leading spectroscopic counterpart to DES. 
It will repeatedly monitor the DES SNe field s using the 2dF 
mult i-object fibre spectrograph (AAOmega, ISaunders et alj 
l2004f) at the Anglo-Australian Telescope (AAT). Its main 
science goal is to measure the redshifts of Type la SNe host 
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galaxies. In addition, a number of fibres in each field will be 
dedicated to monitoring a selected group of quasars to per¬ 
form reverberation mapping. OzDES, in conjunction with 
DES, will monitor ~ 500 quasars for the full five years of 
the survey over the redshift range 0 < z < 4. This is ap¬ 
proximately a ten-fold increase in number and redshift range 
over the existing reverberation mapping samples, and com¬ 
parable in number to the ongoing Sloan Di gital Sk y Surv ey 
(SDSS) reverberation mapping project JShen et alJl2015 4. 

While the OzDES reverberation mapping programs will 
monitor hundreds of quasars, the likelihood that this pro¬ 
gram will successfully recover reverberation lags depends 
on the frequency and accuracy of the light curve measure¬ 
ments, the length of the survey, and the intrinsic variabil¬ 
ity of the monitored AGN sample. OzDES is expected to 
target each field approximately 25 times over the five year 
period. This number of observation epochs is significantly 
smal ler than traditional r everberation mapp i ng campaigns 
(e.g. IPet erson et ~ahl 20021 : iBentz et al.ll2009al : iDennev et al.l 
120091: iBarth et al.l 1201 il l . which have found that emission¬ 
line lag recovery generally requires 30-50 well-spaced epochs 
of observations, and favourable continuum flux variations. 
Our current study aims to investigate the expected efficiency 
of the OzDES reverberation mapping campaign by generat¬ 
ing realistic AGN light curves and attempting to recover 
the input lags. We will then use our findings to determine 
how to optimally select our target AGN sample, and make 
predictions about the scientific results for our sample. We 
also investigate ways to improve the program design and 
execution, through increased cadence, changes to the survey 
length, and improved measurement accuracy, as a means to 
maximise the scientific output. 

The outline of the Chapter is as follows, a technical sum¬ 
mary of the DES and OzDES surveys is given in Section [21 
followed by a description of the survey simulation in Section 
[3] The predictions of the efficiency of the survey are pre¬ 
sented in Section [I] along with the expected improvements 
for several possible survey extensions. We examine the pre¬ 
dicted scientific results in Section [5] Finally, the results are 
summarised in Section [S] 


2 DES/OZDES SURVEY 
2.1 Fields 

The DES SNe fields were chosen to have extensive past ob¬ 
servation histories and to overlap with the Visible and In¬ 
frared Survey Telescope for Astronomy (VISTA). The ten 
chosen fields were the Elias fields (El, E2), SDSS Stripe 82 
(SI, S2), the Chandra Deep Fields (Cl, C2, C3) and the 
XMM Large Structure Survey fields (XI, X2, X3). The co¬ 
ordinates of the fields are given in Table [1] 


Table 1. DES SNe fields 


Target name 

RA (h m s) 

Decl. (° ”’) 

Type 

El 

00:31:29.9 

-43:00:34.6 

shallow 

E2 

00:38:00.0 

-43:59:52.8 

shallow 

SI 

02:51:16.8 

+00:00:00.0 

shallow 

S2 

02:44:46.7 

-00:59:18.2 

shallow 

Cl 

03:37:05.8 

-27:06:41.8 

shallow 

C2 

03:37:05.8 

-29:05:18.2 

shallow 

C3 

03:30:35.6 

-28:06:00.0 

deep 

XI 

02:17:54.2 

-04:55:46.2 

shallow 

X2 

02:22:39.5 

-06:24:43.6 

shallow 

X3 

02:25:48.0 

-04:36:00.0 

deep 


Note: RA and Decl. are given for J2000. 


the KX methocQ (IWarren. Hewett fe Folt j|2000l j using data 
from DES and the VISTA Hemisphere Survey (VHS, 
iMcMahon et al.ll2013l : [Sutherland et aljfeoisl ). and; 3) point 
sources with m riPa f <21 mag selected through photo-z tem¬ 
plate fitting using DES, VHS and Wid e-field Infrared Survey 
Explorer (WISE. [Wright et al.l l201oh photo metry. The de¬ 
tails o f this selection process are described in lBanerii et alJ 
J2015ll . We obtained spectra of 3331 quasar candidates in 
2013, and after visual inspection, the sample was reduced to 
the current 989 objects. Fig.[T]shows the redshift and r-band 
magnitude distribution of this sample. We have ranked the 
sources based on the quality of the spectra and the number 
of emission lines present in the spectra. The highest priority 
objects are shown in blue. When we decrease this sample to 
50 per field, we will incorporate the results of the simulations 
presented here into our target selection criteria and also se¬ 
lect for the most variable quasars (for which lag recovery is 
most probable) based on the first two years of data. 


2.3 Photometry 

DES uses the DECam instrument (Flaugher et al.l 2010 1 
on the Blanco 4-meter telescope at the Cerro Tololo Inter- 
American Observatory in Chile. Images of the ten SNe fields 
are planned to be taken approximately every 5—7 days be¬ 
tween September 1 and February 15 every year between 2013 
and 2017. Photometry from November 2012 - February 2013 
were also taken as part of the DES Science Verification pe¬ 
riod. For a given field, images are taken in all four filters 
in the same night, when possible^. Otherwise, images in the 
remaining filters are taken during the next available night. 
The approximate observation period of DES is shown in Fig. 
[2] along with the visibility of the 10 DES SNe fields. Dur¬ 
ing every night of DES imaging, the SNe fields that have 
not been observed in the last five nights are given highest 
priority, with special preference for the deep fields, C3 and 
X3. To date, the median gap between consecutive observa¬ 
tions is 6.5 days and the maximum gap ranges between 12 to 
21 days. For these simulations we assume that photometric 
observations are taken every 7 days. The nominal exposure 


2.2 Target selection 

The quasar candidates were initially chosen from: 1) known 
quasars in the DES SNe fields with m rjPa f < 21.2 mag; 
2) point sources with m IiPS f < 21 mag selected through 


1 This method selects quasars based on excess flux in the K band 
relative to stars, which is due to the power law nature of the 
quasar spectral energy distribution. 

2 All four filters are observed in one night more than 80% of the 
time 
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Figure 1. The redshift (left) and magnitude (right) distributions of the current OzDES quasar sample. The priority sample consists of 
the highest quality candidates based on visual inspection of the spectra. 


Table 2. DES exposure times and limiting magnitudes of the SN 
fields. 


Shallow Field Deep Field 

Filter Exposure Limiting Exposure Limiting 
Time (s) Mag (AB) Time (s) Mag (AB) 


g 

175 

24.9 

600 

25.6 

r 

150 

24.3 

1200 

25.4 

i 

200 

23.9 

1800 

25.1 

z 

400 

23.8 

3630 

24.8 


times and corresponding limiting magnitudes for both the 
deep and shallow fields are given in Table [7] 

2.4 Spectroscopy 

The OzDES spectroscopic observations are being taken with 
the AAOmega spectrograph fed by the Two Degree Field 
(2dF) multi-object system on the AAT. The 2.1 degree di¬ 
ameter field of view of 2dF is almost identical to that of DE- 
Cam, making it the ide al instrument for the spectroscopic 
followup of DES targets llYuan et al.ll2015l) . and OzDES will 
run over a similar period of time to DES. The 2dF multi¬ 
object system is a robotic fibre positioner that allows simul¬ 
taneous observations of up to 392 targets anywhere within it 
field of view. The projected fibre diameter of the instrument 
is approximately two arcsec. The 2dF fibres feed AAOmega, 
a double beam spectrograph with a wavelength coverage of 
3750A-8900A and a resolution of R ~ 1500. 

During each run, OzDES integrations will be 2 hours 
long in each DES field. Over the five years of the survey, a 
total of 100 nights have been allocated for OzDES with a 
graduated allocation plan. A larger number of nights have 
been allocated each year as the survey progresses: 12 nights 
in 2013, 16 in 2014, 20 in 2015, 24 in 2016, and 28 in 2017. 
It takes approximately four nights to observe all 10 DES SN 
fields (including a 33% allowance for bad weather). We ex¬ 


pect that each field will be visited approximately 25 times 
over the five year period, and the cadence of the measure¬ 
ments will be approximately monthly within each year. The 
2014 observations are now complete and each field has five 
or more epochs of measurements. To date, we have been able 
to devote fibres to 100 AGN per field to help optimise our 
final sample selection. This fibre allocation is scheduled to 
be reduced to a final 50 per field in 2015. 

To give insight into the quality of the current OzDES 
spectra, the distribution of preliminary line flux signal-to- 
noise ratio (SNR) measurements is shown in Figure [3] The 
median SNR value for the total sample is SNR a n ~ 11 and 
for the priority sample is SNR pr i or it y ~ 21. More details 
on the calculation of these preliminary SNR values and ex¬ 
amples of reduced un-flux calibrated OzDES spectra across 
the redshift and magnitude range of the survey are given 
in Chapter ??. OzDES is also monitoring 10-15 F stars per 
field. These observations are important for accurate abso¬ 
lute spectrophotometric calibration of the spectra, which is 
expected to be g ood to <10% based on the results from the 
GAMA project dHopkins et al.|[2013l ). 


3 SIMULATION SETUP 

Our ability to accurately recover lags is highly dependent 
on the presence of prominent features in the continuum and 
emission-line light curves. Accurate detection and character¬ 
isation of light curve features depends on (1) the frequency 
of measurements, (2) the length of the survey, (3) the accu¬ 
racy of the measurements, and (4) the intrinsic variability of 
the object. In order to predict how well OzDES will be able 
to recover lags, we analysed mock catalogues of quasars with 
realistic continuum and emission line light curves. Below we 


4 http: //ww. dartmouth. edu/~physics/labs/skycalc/f Iyer. html 
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Figure 2. The number of hours the DES supernova fields are 
visible throughout the year with an airmass of < 2. The values are 
calculated for the period April 2014 - April 2015 using JSkyCal^. 
The shaded observation period roughly represents the time when 
photometric and spectroscopic data will be taken with the current 
DES and OzDES program design. Section 14.21 investigates the 
improvements afforded by extending this observation peri od to 

fully encompass the time when the fields are visible. 

Priority ■ ■ ■ ■ 



0 20 40 60 80 100 

SNR 

Figure 3. The distribution of emission line signal-to-noise ratios 
(SNR) in the current OzDES quasar sample for the full and prior¬ 
ity sample (defined in section T2.2H . The emission line SNR values 
were estimated for H/3, MgII, and Civ, in each object for all 
observation epochs, using the reduced un-flux calibrated OzDES 
spectra. However, the SNR assigned to each object corresponds 
to the median SNR value of the emission line with the highest 
median SNR value. 


describe our methods for simulating and subsequently re¬ 
covering lags from such light curves. 


3.1 Mock Catalogue 

We began by constructing mock catalogs of 520 AGN dis¬ 
tributed uniformly in 40 redshift bins and 13 magnitude bins 
over the range 0 < z < 4.0 and 18.0 < r < 20.5. This range 
roughly corresponds to the capabilities of the 2dF spectro¬ 
graph within the framework of the OzDES program design. 
The small number of AGN brighter than 18th magnitude 
will be automatically targeted and are not considered for 
this investigation. 


3.2 Monochromatic continuum luminosity 
estimation 


The observed R — L relationship for each broad emission 
line is constructed empirically using the monochro¬ 
matic luminosity of a nearby continuum region as a 
proxy for the ionising luminosity. The most common 
emission-line-continuum region pairs used for con- 
structi ng R — L relat io nship s, are H/3 with L^ lnn ^ 
(e.g. IPeterson fe Hornel 12004 : _ IVestergaard fe Petersonl 


120061: Bentz et al. 2009al. 2013|), Mgn with L^mn A 

( Vestergaard fe Osmerll2pQ9l; Trakbtenbrot fe Netzeil 201211 . 


and C iv with L-i^a (IVestergaard 2002 : IPeterson et al.l 


20041: IV estergaard fc Pe tersonl 20061: iKaspi et all 2007: 


Vestergaard & Osmer 20091 : Park et al.ll2013i : iTrevese et alj 

20141) . and we used these pairs in our analysis. As we are 
restricted to the observed optical spectrum, the redshift 
range is broken into five different sections according to the 
emission lines we are able to observe: z < 0.54 H/3 only, 
0.54 < z < 0.62 H/3 and Mg II, 0.62 < z < 1.78 Mgn 
only, 1.78 < z < 1.96 Mgn and Civ, and 2 > 1.9 Civ 
only (see Fig. [1] for illustration). Accurate measurement 
of the surrounding continuum is required to measure the 
line flux correctly. Therefore, these redshift ranges allow a 
generous amount of continuum on either side of the broad 
line, and avoid bluer wavelengths where the spectrograph’s 
throughput is low. However, our choice of redshift ranges 
are somewhat conservative with respect to how far into the 
blue and/or red each emission line or continuum region can 
be well-measured. 

The monochromatic fluxes at 5100A, 3000A, and 1350A 
for each AGN are estimated using the known DECam r- 
band magnitudes and redshifts, and a K-correction based 
on the filter response curves of the DECam r-band filter 
and the Sloan Digital Sk y Survey (SDSS) quasar template 
dVanden Berk et al J200ll) . The SDSS quasar template spans 
the rest-frame wavelength range 800A-8555A and is based on 
a sample of quasars that covers a similar redshift and mag¬ 
nitude range to that being observed by OzDES, so it should 
be a representative template for this work. The quasar tem¬ 
plate is scaled to the input magnitude of each target un¬ 
der the assumption tha t the bolometric luminosity scales as 
Lt,oi = 9.0ALa( 5100A) (iKaspi et al.ll2000l ). This bolometric 
co rrection is the approxi mate midpoint betwee n the r esults 
of iRichards et al.l ll2006l) and iKrawczvk et al.l ll2013l ) who 
found values of 10.3 ±2.1 and 7.79 ± 1.69, respectively. The 
luminosity is then estimated assuming a ACDM cosmology 
with Ho = 70 km s -1 Mpc -1 , Qm = 0.3, and 12 a = 0.7. 
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2000. 4000. 6000. 8000. 10000. 

Wavelength(A) 

Figure 4. The five redshift regimes in which we are sensitive to 
different broad lines (or pairs of lines) given the properties of the 
2dF instrument. The black curves show the observed frame wave¬ 
length of each broad emission line considered in our simulations 
as a function of redshift. 


This process does not take into account host galaxy 
contamination, Galactic extinction, or intrinsic variations in 
the AGN spectral energy distribution (SED). Nonetheless, 
it provides a sufficient luminosity estimate for our present 
purposes. When host galaxy contamination is included in 
the simulations, no significant change in the results was ob¬ 
served because the loss in sensitivity due to contamination 
from the host galaxy is generally balanced by the increased 
variability expected for fainter objects (see section l3~4l) . For 
more discussion on this point see Section [5] 


3.3 Lag estimation 


We calculate the time lags associated with each source from 
the monochromatic continuum luminosity and previously 
published R — L relationships for H/3, Mg II and CIV. We 
adopt the H/3 R — L relationship 


log 10 (Rh? It days) = - 21.2 ± 2.2 + 0.517 ± 0.033 (2) 

x log 10 (AL a ( 5100A) erg s -1 ), 

presented bv lBentz et al.l (l2009al ). which is derived from RM 
measurements of 35 AGN spann ing four orde r s of m agnitude 
in luminosity. Switching to the lBentz et al.l 120131) relation 
would not appreciably change our results. The CIV R — 
L relationship is not as well determined as the H/3 R — L 
relationship b ecause CIV lags have only been measured for 
a few objects ( Koratkar & Gaskclll fl99ll : IPeterson fe Hornel 
120041 ; iKaspi et aL 20071 ). Most of these objects have similar 
luminosities, and there is only a single object at each of the 
low and high luminosity ends. The CIV R — L relationship 
we use, 


log 10 (i?crv It days) = - 23.3 ± 2.6 + (0.55 ± 0.04) (3) 

x log 10 (ALa( 1350A) erg s -1 ), 


is taken from lKaspi et all J2007I ) and is based on only seven 
objects. Mg 11 is not yet well-studied with reverberation 
mapping, and there are only a few studies in which the time 
lag has been measur ed jMetzroth, Onken fe Pet ersonl 120061 : 

I Reichert et al.l Il994l : iDietrich fe Kollatschnvl 1 19951 . though 
the latter two papers only present marginal detections). 
However, because the Mgn lag measurement was found 
to be consistent with t h e H/3 lag in two of the stud¬ 
ies jReichert et al.l 1 19941: IStirpe et al.l 1 19941; iBentz et alj 


I2006L using iMetzroth. Onken fe Petersonl ( 2009 )). there is 
a strong correlation between the width of Mgn and H/3 
emission lineal and H/3 and Mgn have similar ionisation 
parameters, it is generally assumed that the two lines 
originate at the same radiu s from the ionising source. 
iTrakhtenbrot fe Netzeij ((20121 ) used this assumption to esti¬ 
mate a Mg 11 R — L relationship of 


log 10 (R Mg ii It days) = - 25.72 ± 0.62 + (0.615 ± 0.014) 

(4) 

x log 10 (AL a ( 3000A) erg s -1 ). 

This estimate was created using an empirical correlation 
between ALa( 5100A) and ALa( 3000A), and an existing H/3 
R — L relationship. As a consequence, it is less certain than 
an R—L relationship derived from direct reverberation map¬ 
ping measurements, but it is our only option in the absence 
of such measurements for Mgn. 

The resulting distributions of the observed time lags 
(including time dilation) with redshift and magnitude are 
shown in Fig. [S] The Civ lags tend to be a factor of ~ 2 
smaller than the Mgn lags f or the same object s. This is 
consistent with the findings of lKaspi et alj d2007l ). 


3.4 Light Curve estimation 

We model the AGN continuum light curves as a damped 
random walk (DRW) characterised by a time scale, 7X0, 
and the asymptotic ampli tude of th e struc ture function over 
long time scales, Sp(oo). IZu et all (I 2 OI 3 I ) show this model 
is a good representation of quasar variability on the time 
scales of the OzDES survey. Our ability to recover the time 
lag for each AGN is highly dependent on these variability 
parameters. If Sf(oo) is large and td is short, it is more 
likely that significant variations in the light curve will be 
seen du ring the observatio n period. 

iMacLeod et al.l (l2010t) found that the values of td and 
Sf(o o) scales with the luminosity of the AGN, the observed 
wavelength, and the mass of the black hole, following the 
power law, 

logioC®) + B a log 10 (A/4000.) + C a {Mi + 23) 

+ D a log 10 (M BH ), (5) 

where a refers to either td (in days) or S'f(oo) (in mag), 
A (in A) is the continuum wavelength of interest (1350A, 

5 However, this correlation is not a tight one-to-one relationshi p 
between the Mgn and H/3 width dTrakhtenbrot &; Netzej|2012l) , 
as is commonly reported in the li teratu re. 

6 Note that IMacLeod et alj (|2010|) use r without the 
subscript for this parameter, but we use td following 
IZu, Kochanek &; Petersonl <l201lh so as not to confuse this pa¬ 
rameter with the reverberation lag, which we denote as r. 
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m r z 

Figure 5. The distribution of lags in a mock AGN sample, as a function of r-band magnitude (left) and redshift, z (right). The different 
colours and symbols correspond to the different emission lines (black circles: H/3, red triangles: Mgn, blue crosses: Civ). 


3000A or 5100A), the coefficients are 


A td = 2.4, 


B td = 0.17, 


C TD = 0.03, 


D td = 0.21, 


A Sf = -0.51, B Sf = -0.48, C Sf = 0.13, D Sp = 0.18, 

Mi is the absolute magnitude of the quasar and is calculated 
using the kn own r-band ma g nitud e and the K-corrections 
prescribed in iRichards et all d200(t for SDSS, and Mbh is 
the mass of the black hole. The black hole masses were ran¬ 
domly assigned using the probability distribution described 
by IMacLeod et al.l (1 20101 1 , 


P(log 10 = 


■ exp 


(logic Mbh — log 10 Mbh ) 2 


2 o m bh 


where log 10 M 


y/2ncr 

(6) 

bh = 2.0 — 0.27 Mi and the spread 

in the distribution, &m bh = 0.58 + 0.011 Mi. The 

mass of the black hole is in Solar units. The DES 
and SDSS qri filte rs are quite similar, so the use of 
the Richards et al.| (f2 006ll K-corrections are j ustifie d 
llHonscheid. DeP ov fe for the D ES Collaboration! l2008li . 
Additionally, MacLeod et al.l (j2010h assumed the same cos¬ 
mological model as we use here, so the absolute magnitude 
scales are identical. 


3.4-1 Continuum light curve 

The continuum light curve, in magnitude, is defined as a 
combination of a mean g. and variable term A C(ti), such 
that G(f) = fx + A C(t). The value of g. is defined as the 
monochromatic continuum flux at given wavelength (5100A, 
3000A, and 1350A; Section [3~2ll converted to a magnitude. 
For a DRW, the variable component is constructed by ini¬ 
tialising the light curve at to as AG(to) = aG( 1) where 
G(l) is a Gaussian random deviate of unit dispersion, and 
Sf( 00 ) = \/2a. Sub sequent points are created using the re- 


cursion formula (see 

Kcllv, Bechtold & Siemieinowska 20091: 

IMacLeod et al.ll201C 

: Kozlowski et al. 2oTo|) . 


AC(t i+ i) =A C(U) exp (-|t;+i - U\/td) (7) 

+cr [1 - exp (—2|t i+1 - ti\/r D )] 1/2 G( 1). 


3-4-2 Emission line light curve 


The emission line light curve is the response to the contin¬ 
uum light curve. If we only consider the temporal response 
of the overall emission line, the resulting light curve is given 
by 

A L{t)= [ $(r)A C(t-T)dr, (8) 


where A L(t) is the emission line light curve flux rela¬ 
tive to its mean value, AG(f) is the variable component 
of the continuum light curve defined above, and T(t) is 
the transfer function, which describes the emission line 
response to a delta function outburst in the continuum 
jBlandford fe McKeelll982h . The transfer function is related 
to the overall structure and kinematics of the broad line re¬ 
gion and its true form remains an active area of research 


(Pancoast. Brewer & Treu 

l201ll; Zu. Kochanek & Petersonl 

I201II: ’Pancoast et al. 2012 

: iGrier et al.ll2013lb For simplic- 


ity we hav e chosen a top hat transfer fun ction, mirroring 
JAVELIN (IZu, Kochanek fe Peterso~71l201llh given by 


T(f) 


A 

2w 


0 


for T — W < t < T + w 
otherwise 


(9) 


where A is a scaling term and w is the half-width of 
the top hat function. We set the scaling term, A = 1, 
and the half-width to ui = O.lr. This half-width value 
was m otivated by previ ous reverberation mapping cam¬ 
paigns dGrier et al.ll20l3 l. The choice of amplitude makes 
a strong assumption about the line’s responsivity to con¬ 
tinuum variations. Photoionisation research suggest that 
diffe rent lines will have different assoc i ated responsivi- 
ties (iGoad. O’Brien fe Gondhalekail 1 19931 : iKorista fc Goadl 
l200d . 2004), and there is observational evidence that the 
responsivity for individual lines is not consistent be tween 
monitoring program s (e.g. Mgn responsivity: IWool 120081 : 
ICackett et al.|[2015h . Additionally, the value of A C(t) imple¬ 
mented in this analysis is only a proxy of the true ionising 
flux and as a consequence, the transfer function quoted will 
be different to the inherent transfer function of the system. 
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We investigate the effects of our choice of transfer function 
in Section [6] 

3-4-3 Measurement uncertainty 

We created mock light curves with daily cadence, starting 
sufficiently before the start of the observing campaign to 
allow time for the emission line response (described by the 
convolution in Equation ®) to lie in the observed time pe¬ 
riod. We then down-sample these light curves following the 
OzDES program design and any extensions described below. 

A Gaussian error is added to each mock light curve 
measurement based on the expected flux uncertainties. The 
expected flux measurement uncertainty for OzDES R.M 
project is 0.01 mag for photometry and 0.1 mag for spec¬ 
troscopy. We do not include magnitude dependent (photon 
counting) errors on the line flux measurements as our un¬ 
certainties are generally expected to be dominated by our 
overall absolute flux calibration. The absolute flux calibra¬ 
tion of the spectroscopy is performed using the following 
procedure: 1) the observed F stars are matched with a F 
star with equivalent g-r colours (measured from DES pho¬ 
tometry) from an existing stellar catalogue; 2) the catalogue 
spectrum is then warped to exactly replicate the colours of 
the observed F star; 3) next, the observed spectrum is di¬ 
vided by the warped catalogue spectrum; 4) the resulting 
function is smoothed and represents the sensitivity curve 
of the observed spectrum. This process is performed sep¬ 
arately for the red and blue arm; 5) and repeated for all 
the observed F-stars; 6) the median sensitivity is then com¬ 
puted for each season. The sensitivity is dependent on the 
observed wavelength and radial position of the fibre, and 
the median scatter in the sensitivity curve is ~5%; 7) At 
that time, all AGN spectrum are corrected according to the 
measured sensitivity curve and synthetic g , r, and i magni¬ 
tudes are calculated; 8) Finally, the spectrum is compared 
to the nearest photometric g, r, and i magnitudes (almost 
all within four days) and scaled under the assumption that 
there will be insignificant changes in flux on this time scale. 
The photometric uncertainty in this case is on order of a few 
percent. Current OzDES data shows that residuals of these 
fits are roughly Gaussian in form, therefore it is reasonable 
to assume Gaussian uncertainties in our light curve measure¬ 
ments. For realism, we further randomly shifted 10% of the 
data points without adjusting their error bars to introduce 
an element of the non-Gaussianities present in real data. 
The shift was chosen from a Gaussian distribution with a 
mean value of zero and a standard deviation twice the size 
of the measurement uncertainty. An example of a simulated 
light curve is shown in Fig. [6] 

The accuracy of the emission line flux measurement also 
sensitive to the method of determining the emission line 
flux. The method for measuring the emission line flux in 
the OzDES data will differ depending on the quality of the 
spectrum. If the signal-to-noise ratio in the individual spec¬ 
trum is sufficient, the spectrum will be modelled to separate 
the emission lines of interest from the continuum, absorp¬ 
tion line and contaminating emission line signal. Otherwise, 
the emission line flux is simply calculated as the integrated 
flux above an approximate linear continuum fitted between 
two pseudo-continuum regions near the line. Mgn A2798 
flux measurements are especially affected by Fe II contami¬ 


nation (IVestergaard fo WilkestfioQlh and therefore care must 
be taken when calculating the flux for this line. Additionally, 
AGN with strong absorption (e.g. broad absorption lines; 
BALs) will be discarded from the final sample. 


3.5 Possible Survey Extensions 

We also consider various extensions to the baseline survey 
to determine if modifications to the program design can in¬ 
crease the scientific return. These extensions include (i) de¬ 
creasing the effective seasonal gap; (ii) increasing the ca¬ 
dence; (iii) improving the data quality; and (iv) extending 
the total survey duration. While ultimately uncontrollable, 
we also consider the effect of a reduction in the number of 
observed epochs due to extreme weather losses. 

Seasonal Gap [Full season/Year] : We consider two 
extended spectroscopic observation windows to investigate 
the effect of the seasonal gap on lag recovery. The longest 
possible observation window in which any one field can con¬ 
tinuously be observed is approximately between May 1 and 
Feb 14 for fields El and E2 (see Fig. [2j. This will be our first 
possible extension and will be referred to as ‘full season’. As 
the first two years of observations are completed or under¬ 
way, we can only apply this extension to the following three 
years. We also consider a case in which observations are 
taken over the full year, with no seasonal gaps. This is our 
second possible extension and will be referred as ‘year’. Due 
to the restricted observability of the DES SN/OzDES fields 
resulting from their position in the sky, this observational 
setup is not possible for OzDES, but illustrates the poten¬ 
tial for future surveys targeting polar fields. In both cases, it 
is assumed that photometric measurements are taken weekly 
and spectroscopic measurements are taken monthly. 

Cadence [Weekly]: We also consider the case where 
spectroscopic measurements are taken weekly over the base¬ 
line observation window of the DES photometry as defined 
in Section [2] This will test how the spectroscopic measure¬ 
ment frequency affects the recoverability of the lag. Previ¬ 
ous reverberation mapping campaigns have found that high 
quality and high cadence in the continuum light curve is 
required for accurate lag recovery as it is the driver of the 
line light curve. However, it is also important that the line 
flux measurements are taken with sufficient cadence to map 
the line response accurately. This possible extension will be 
referred to as ‘weekly’. Again, we can only apply this exten¬ 
sion to the last three years of the survey. We should note 
that this extension increases the number of spectroscopic 
measurements signifi cantly, which also plays a major role in 
the recovery of lags dHorne et al . 20041 : 1 Shen et"afll2015h . 

Data quality [Goal]: We investigate the effects of data 
quality by changing the spectroscopic measurement uncer¬ 
tainty. We use the baseline sampling rate of OzDES and 
reduce the spectroscopic measurement uncertainties from 
0.1 mag to 0.03 mag, which corresponds to the optimistic 
goal for the calibration of the OzDES spectral data. The 
photometric measurement uncertainty is kept constant at 
0.01 mag. This extension will be referred to as ‘goal’. We 
also consider a spectroscopic measurement uncertainty of 
0.03 mag combined with the ‘year’ extension, referred to as 
‘year+goal’. 
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Survey Length [Long]: Finally we investigate an ex¬ 
tension of the survey by 2 years with sampling and cadence 
equivalent to the last year of the planned DES/OzDES sur¬ 
vey (i.e. 7 spectroscopic measurements per additional year). 
This possible extension will be referred to as ‘long’ and it 
should allow the recovery of longer lags and more accurate 
recovery of lags for a broader AGN population. 

Weather: The expected 25 epochs, used in the simula¬ 
tions, already takes into account the expected weather loss 
for the A AT. However, to consider the effects of extreme 
weather we simulated losing an additional 3-5 spectroscopic 
epochs over the 5 year period. 


150 iterations per chain for both the burn in periocQ and to 
sample the posterior probability distribution. This is equiv¬ 
alent to 22500 total burn in iterations. For the case where 
both emission lines light curves are fitted simultaneously, we 
used 200 MCMC chains with 200 iterations per chain. The 
larger number of iterations accounts for the additional com¬ 
plexity in the double line fit. For the baseline OzDES setup 
we analysed ten realisations of the mock catalogue of 520 
AGN. For the investigation of different survey extensions, 
we only analysed two realisations of the mock catalogue per 
extension. 


3.6 Recovering time lag 


Traditionally, reverberation lags have been recovered us¬ 
ing sim ple linear interpolation and cross correlation tech¬ 
niques jjGaskeUfcS park^ Il986l : IWhite fe Petersonl 1 1994 
lEdelson fc Krolikl 1988l l. but recently other approaches have 
been implemented that take into account our existing knowl¬ 
edge of AGN behaviour as a means to optimize, and in some 
cases, to improve the likelihood of accurate lag recovery. We 
adopt one such approach, using the program JAVELIN, an 
upda ted version of SPEAR (see IZu. Kochanek fc Petersonl 
l201ll. for details). Instead of linearly interpolating between 
data points, JAVELIN uses a damped random walk to model 
the AGN continuum light curve, and attempts to fit the 
emission line light curve by convolving the continuum light 
curve with a top hat transfer function. 

JAVELI N uses the amoeba minimisation method 
(IPress et al.l 1 19921 1 to recover a model of the contin¬ 
uum light curve (including DRW parameters) and trans¬ 
fer function that best fits the continuum and the line 
data. JAVELIN has been generally found to be con- 
sisten t with traditional cross- c orrela ti on (CCF) meth ¬ 
ods dZu . Koch anek fc Petersonl 1201 ll : iGrier et~al ] 120121: 
IPeterson et al.ll2014l 'l and has the advantage of fitting mul¬ 
tiple emission line light curves at once. The biggest limi¬ 
tation of JAVELIN is that its error estimates assume well- 
characterised Gaussian noise, so caution must be used when 
interpreting the parameter uncertainties if these assump¬ 
tions are violated. 


3.6.1 Implementation 

We allow JAVELIN to explore the lag range between zero 
days and three times the input lag, with a maximum allowed 
lag of 1931 days, corresponding to the separation between 
the first photometric and last spectroscopic measurement 
planned for OzDES. Lags longer than 1931 days cannot be 
constrained by the data; however, we expect very few, if 
any lags this large in the considered magnitude and redshift 
range (see Figure O. The window of 0 < t < 3r was found 
to be sufficient to fully enclose the recovered likelihood dis¬ 
tribution of the lags, and there is no evidence for an artificial 
cutoff in the likelihood distribution (see Fig. [6]). 

We analysed the five different redshift ranges separately. 
For the redshift ranges where two emission lines are present, 
the two lags were fit both individually and simultaneously. 
For the single line case we used 150 MCMC chains with 


3.6.2 Output 

For each continuum and emission-line light curve, JAVELIN 
produces an estimate of the posterior distribution for all the 
fitted parameters, including the lag. In a substantial number 
of cases, a clear peak is present and easily identified in the 
lag posterior distribution, corresponding to the best fit lag 
value. In the remaining cases, multiple peaks of comparable 
size were present in the distribution or in some rare cases 
no distinct peak was detected at all. 

We devised a simple method to classify each posterior, 
trained on a manually classified sample of posteriors. In 
the most basic terms, we smooth the posterior distribution, 
and identify the dominant peak and any secondary peaks. 
We classify the quality of the lag based on the existence 
and relative size of the secondary peaks, which corresponds 
to the relative probability of the lag being associated with 
each peak. We classified the output probability distributions 
into four categories: Accepted - grade 1, Accepted - grade 
2, uncertain, and rejected. Accepted - grade 1 means that 
no secondary peaks were present in the posterior distribu¬ 
tion. Accepted - grade 2 means that a secondary feature 
is present but the ratio between the probability of the pri¬ 
mary peak and relative probability of the secondary peak 
is smaller than 15%. The secondary peak may be a bump 
associated with the main peak. The relative probability of 
the secondary peak is defined as the difference between the 
maximum probability value of the secondary peak and ei¬ 
ther zero or where the secondary connects to the main peak, 
whichever is the smallest value. This definition is chosen 
to allow small deviations in the probability distribution of 
the main peak. Uncertain posteriors had multiple peaks of 
similar magnitude, while rejected posteriors did not have 
any clear maxima. The classifications were made without 
prior knowledge of the true lag. Some examples of typical 
accepted, uncertain, and rejected lag posterior distributions 
are shown in Fig. [6] Once classification has been performed, 
we redefine the lag window boundaries for the Accepted - 
grade 2 cases to exclude the secondary peak where applicable 
and calculate the credible regions around the primary peak 
only. This algorithm performed reasonably well compared 
to human classification, with a 5% misidentification rate. In 
practice, the real data and lag fits will be subjected to man¬ 
ual inspection which will reduce the number of misidentified 


7 The burn in period refers to an initial portion of a Markov chain 
sample that is discarded to minimise the effect of initial values 
on the posterior inference. 
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lags. The recovered lag used for the rest of this Chapter is de¬ 
fined as the median lag value from the restricted JAVELIN 
Monte Carlo chain for accepted lags only, and the lag uncer¬ 
tainty is given by the 16th and 84th percentile range about 
the primary peak. 


3.6.3 Performance metrics 

To quantify the performance of the OzDES reverberation 
mapping campaign and its possible extensions, we define 
three performance metrics: 1) the recovered fraction, describ¬ 
ing how many lags we recover, 2) a a, quantifying how accu¬ 
rately the recovered lags are measured, and 3) the number 
of misidentified lags in the recovered fraction. 


Recovered Fraction: Using the classification described 
in the previous section for accepted, uncertain and rejected 
posteriors, we define the recovered fraction as the fraction of 
Accepted-grade 1 and Accepted-grade 2 lags. We determine 
the natural spread in the recovered fraction by performing 
bootstrap resampling and calculating the recovered fraction 
for 1000 iterations of the sample. The stated recovered frac¬ 
tion is given as the median recovered fraction from the re¬ 
sampling and the uncertainty is given by the 68th percentile 
values. 

Accuracy (cta): We next quantify the accuracy of the 
recovered lags. The accuracy is defined by the logarithmic 
ratio between the median recovered lag and the true lag, 

/ rec \ 

A = lo s(^J- (io) 

Fig. m shows that the distribution of A values is approxi¬ 
mately Gaussian, for both the single and two-emission-line 
simulations, and the skewness of the distributions is mini¬ 
mal. The underlying width of the cta distribution charac¬ 
terises the accuracy with which we can expect to recover 
lags. As each recovered lag has its own uncertainty, we can¬ 
not simply calculate the unweighted standard deviation. In¬ 
stead, we find the weighted standard deviation by finding 
the value of cta that minimises the likelihood, 


£(A,et|crA) = Y 


Li VM°A+e?) 


exp 


—A? 


L2(^|+e?) 


( 11 ) 


where = {t^ % — is the average measurement 

uncertainty associated with each lag measurement. To find 
the uncertainty on cta we again use bootstrap resampling. 

Misidentified lags: A lag is designated as misidentified 
if it is 3o> away from the true lag, where ay is the estimated 
uncertainty in the lag. This is equivalent to a for symmetric 
uncertainties around the true lag. However, as the uncertain¬ 
ties are asymmetric, we consider the error estimates towards 
the true lag value. Note that while we count the number of 
misidentified lags, we do not remove them from our sample 
when calculating cta or any further analysis. This is because 
we will not know the true lag length when we analyse the 
real data, and we expect misidentified lags will contaminate 
the true sample at some level. 


3.6.4 Two line fitting - Simultaneous vs. Individual 
Fitting 

As previously mentioned, when two emission lines are 
present in the spectrum we fit the corresponding lags both 
individually and simultaneously. Simultaneous fitting has 
the potential to constrain the continuum light curve more 
tightly, which could enable a better recovery of the two lags. 
However, we find the accuracy of the recovered lags and the 
fraction of misidentified lags to be consistent between the 
two methods, although the simultaneous fits had a lower ac¬ 
ceptance rate than the individual fits (Fig. 0. In general, 
the lag posterior distribution for the simultaneous fitting 
case showed a peak at the same location as the individ¬ 
ual fits, but other secondary peaks were sometimes present, 
making the lag classification more uncertain. It is likely that 
this occurs when the lag signal is weak, leading a spurious 
peak in one lag to amplify an otherwise weak peak in the 
other lag posterior, leading to multiple peaks. The reverse 
also occurred, although it was less common, and in this case, 
the information from the simultaneous fitting nullified the 
spurious signal. We use the results from the individual line 
fits for the rest of the analysis. 


4 RESULTS 
4.1 Basic Setup 

4-1.1 Recovery Fraction 

We first investigate the recovered fraction in bins of red- 
shift and magnitude. This will help to optimise our target 
selection choices. The resulting fractions of accepted lags are 
shown in Figures [8] and [9] There is a clear trend in favour 
of low redshift, faint objects. 

Since the observer-frame lag also depends on magni¬ 
tude (through the R-L relationship) and on redshift (due to 
time dilation), we also investigated how the recovered frac¬ 
tion correlates with the input observed frame lag. Fig. [TD] 
shows the resulting fraction of each lag quality classification 
(accepted, uncertain, and rejected) as a function of input 
observed frame lag. There is a steady decline in the recov¬ 
ered fraction with lag length, and a reciprocal increase in 
the fraction of uncertain and rejected lags. The likelihood 
of successfully recovering the lag increases when the over- 
all survey length is appreciably longer than the lag length 
jHorne et al.ll20(M h because more light curve features can 
be traced by both the continuum and emission line light 
curves. However, program cadence is also important for lag 
recovery, where sampling the light curves at the “Nyquist 
frequency” is required to accurately resolve the lag. Conse¬ 
quently, objects with an expected observer frame lag shorter 
than the minimum emission line light curve sampling rate 
(~ 30 days) are less likely to have a reliable lag estimate. 

The observed dependence of the recovered fraction on 
lag length will also be affected by the magnitude dependence 
of the damped random walk parameters, specifically Sf{ 00 ) ■ 
The Sf{oo) parameter describes the long term amplitude of 
the variable continuum light curve component. If Sf{oo) is 
small, as is true for more luminous objects, it is less likely 
that the observed light curve will vary by a significant de¬ 
gree compared to the noise, and produce a lag measurement. 
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Figure 6. Examples of typical lag posterior distributions (left) for the accepted - grade 1 (top), accepted - grade 2 (middle top), uncertain 
(middle bottom) and rejected (bottom) classifications. The dot-dashed vertical line represents the true lag value. The hatched region in 
the accepted - grade 2 case demonstrates the likelihood region excluded from the final lag estimate. The right panels show the input 
continuum and emission line light curves (circles) and true underlying light curve (red) compared to the weighted average of light curves 
that fit the data well from JAVELIN (black line), and the corresponding dispersion of these light curves (grey shaded region). 
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Figure 7. The relative offset of the recovered lag, A, for the individually fitted (blue circles) and simultaneously fitted (olive triangles) 
samples (left). The error bars represent the inner 68 percentile errors. The distribution of A is shown in the right panel. 
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Figure 8. The recovered fraction in each redshift-magnitude bin for the 10 mock catalogues, as a function of magnitude (left) and 
redshift (right). There is a clear trend both in redshift and magnitude to favour lower redshift, less luminous objects. The dashed vertical 
lines correspond to the bounds of the different redshift ranges for the different emission lines. 


More accurately, the recovered fraction will be affected by 
the mean fractional variation of the continuum light curve, 
Fvar, defined as the rms variab ilit y amp li tude of the contin¬ 
uum magnitude fsee lRodrfguez-Pasc.ual et al.lll997r ). rather 
than the true Sf(oo) value. By definition, Fvar is closely 
related to Sf(oo) as Sf(oo) = V2 a, where a is the long term 


standard deviation in the continuum magnitude and Fvar 
is the fractional standard deviation in continuum flux over 
a finite observation period. The approximate transforma¬ 
tion, 10 °' 4<t ~ 1 + Fvar , can be used when the observation 
period is sufficiently long. This approximation tends to un¬ 
derestimate the input Sf(oo) value by roughly 30% over the 
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Figure 9. The 2D distribution of the recovered lag fraction in 
redshift-magnitude bins. This distribution is similar in shape to 
the distribution of lags (Fig. [5j, indicating that the acceptance 
rate is highly dependent on the intrinsic observed frame lag. 
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Figure 10. The fraction of accepted (blue), uncertain (purple) 
and rejected (red) lags as a function of the observed lag. 


five year observation baseline. AGN with Frms values lower 
than the expected measurement uncertainty of OzDES are 
more likely to have an uncertain lag classification. This is 
illustrated in Fig. [Tl] 

In practice the lag length and light curve variability are 
related through the magnitude dependence of the param¬ 
eters (R-L relationship and Equation [5}, so the recovered 
fraction can generally be characterised simply by the ab¬ 
solute magnitude and intrinsically brighter objects are less 
likely to have a lag measurement. 



1o 8( F VAr) 

Figure 11. The effect of the rms variability amplitude of the 
continuum, Frms j on the fraction of accepted (blue), uncertain 
(purple) and rejected (red) lag likelihoods. In general more vari¬ 
able objects are more likely to have a recoverable lag. 


4-1.2 Accuracy of recovered lags 

Next we look at the accuracy of the recovered lags. Fig. 
EO shows how the accuracy changes with observed frame 
time lag. We divided the sample into accepted-grade 1 and 
accepted-grade 2. The accepted-grade 1 lags were more ac¬ 
curate in general, with a <ta = 0.070 ± 0.006, compared to 
<ta = 0.083 ± 0.004 for the grade 2 sample. However, the 
number of grade 1 lags is less than half the number of grade 
2 cases, so excluding them would have a significant impact 
on our final sample size. For both samples we see a slight 
increase in <ta for r < 30, and r ~ 0.5 years. Lags shorter 
than 30 days are not accurately recovered due to the lim¬ 
ited temporal sampling of the survey. The increase at half 
a year is likely a consequence of the seasonal gaps in the 
survey. Lags can only be accuractely recovered when there 
is some overlap between the observed continuum and emis¬ 
sion line light curve features (e.g. associated rise and drop 
in light curve). Lags close to 365 x n/2 days, where n is 
an odd number, will have very few observed continuum and 
emission line light curve measurements that can trace out 
associated light curve features, leading to larger uncertain¬ 
ties and lower accuracy. Conversely, for lags close to a year, 
associated light curve features are well traced by both the 
observed continuum and emission line light curve, leading to 
a well defined and accurate lag estimate. This effect strongly 
suggests extending our observing season, as we investigate 
in the next section. 


4.2 Extensions 

We next consider how the possible survey extensions affect 
the lag recovery relative to the baseline survey. The main 
results from these various extensions are summarised graph¬ 
ically in Figures [T51I15I where Fig. E] shows the distribution 
of A, Fig. ll4l shows the recovered fraction, and Fig. llSI shows 
a a as a function of the real lag for the different extensions. 




























14 A. L. King et al. 



r reai (Observer Frame) [days] 

Figure 12. The scatter around the true lag, <ta, as a function of 
true observed frame lag for the accepted grade 1 sample (purple 
hashed) and accepted grade 2 sample (black solid). The dotted 
lines correspond to half year intervals, and the dashed lines cor¬ 
respond to full year intervals. There are no accepted grade 1 lags 
recovered for T rea i > 1.5 years. 


Seasonal Gap: Implementing the full season extension is 
expected to increase the recovered fraction to ~50%, which 
corresponds to a 20% relative increase over the baseline sur¬ 
vey. The value of cta also decreases significantly, especially 
for lags at 365 x n/2 days. This is expected, as the smaller 
seasonal gaps reduce both the need for interpolation and the 
range of lags for which common features cannot be traced 
by both the measured continuum and emission line light 
curve. The year extension shows still greater improvements. 
However, the full season and year extensions only show sig- 
nihcant improvement in accuracy for short lags, and no ap¬ 
preciable improvement was seen for long lags. In fact, the 
year long extension accuracy was diminished compared to 
the baseline survey for long lags, as the total duration of the 
baseline survey is extended beyond five years due the inclu¬ 
sion of the DES science verification data. Without extending 
the survey length, long lags remain hard to recover because 
the continuum and emission line light curves still only have 
a few common features. 

Cadence: The recovered fraction and accuracy of the re¬ 
covered lags improved significantly with the weekly exten¬ 
sion. The finer sampling enables superior recovery of shorter 
lags, as can be seen in Fig. 1141 It is expected that some of 
this improvement is due to the increased number of epochs. 
Although, the number of epochs is greater in this scenario 
than the full season extension, and the overall improvement 
in the lag recovery is not as large. This suggests that if ad¬ 
ditional telescope time is awarded in the next three years 
of observations, it should be used to extend the observation 
season rather than to have finer sampling. 

Data quality: Reducing the measurement uncertainties 
relative to the baseline significantly increases the recovered 
fraction and accuracy of the recovered lag. However, im¬ 
proving the measurement uncertainty in the year scenario 
(‘year+goal’), appears to have a negligible effect. In fact, it 


appears as though a reduction of the emission line light curve 
measurement uncertainty to 0.03 mags has a gain equiva¬ 
lent to carrying out the year survey. Therefore, if we can 
reduce the uncertainties on the light curve measurements, 
additional epochs will be of less consequence. Alternatively, 
once we reach the limit of reducing the uncertainties, we can 
improve the scientific outcome with more epochs. 

Survey Length: The long extension shows an increased 
recovered fraction and accuracy primarily for the objects 
with longer lags, as expected. We also see a moderate in¬ 
crease in accuracy for all lags. However, a decrease in recov¬ 
ered fraction and accuracy at half a year is still present due 
to the seasonal gaps. 

Weather: Finally, if 3-5 additional epochs of spectro¬ 
scopic data were lost over the five years, only a slight drop 
in recovered fraction and accuracy is expected. 


5 PROSPECTS 

Maximising the return is not simply a question of maximis¬ 
ing the number of recovered lags. Next we use the results 
of the simulations to optimise target selection for measuring 
black hole masses and constraining the R — L relationship. 
Additionally, we fold in the redshift and magnitude distribu¬ 
tion of the OzDES target quasars into our predicted results. 
Fig. [Tg] shows an example of the recovered lags we expect 
for the baseline OzDES survey if the final 500 AGN were 
chosen randomly from the currently observed 989 candidate 
quasars. The acceptance rate is ~ 35% and the sample cov¬ 
ers a luminosity range of 10 39 < ALa(5100A) < 10 46 with 
mean [i [log [ALa(5100A)] } = 45.0, and extends to redshift 
- 4. 

5.1 Black hole mass measurements 

The uncertainty in a RM black hole mass estimate is, 



assuming the errors on the virial factor, /, the line width, 
AV, and BLR radius, R, are independent. Generally, the 
formal errors reported for Mbh only include the uncertain- 
ties from R and AV and ignore th e uncertainty in / (e.g. 
iBentz et al .1120081 : iGrier et al.ll20l3 l even though it is gen¬ 
erally the largest source of uncertainty in the mass determi¬ 
nation. The typical uncertainty in the reverberation m asses 
due to uncertainty in / is ~ 0.43 dex (IWoo et al.ll2010l ). 

The velocity dispersion, AV, is measured from the emis¬ 
sion line width in either the rms variance spectrum (hence¬ 
forth referred to simply as the rms spectrum) or the mean 
spectrum. The rms spectrum should be used whenever pos¬ 
sible as it is the best representation of the responding (vari¬ 
able) component of the emission, and will therefore provide 
the most accurate mass estimate. However, sometimes the 
rms spectrum is not of high enough quality or the rms pro¬ 
file is unusual in shape, and the in terpretation of wid th is 
not strai ght forward (e.g . NGC4151: IBentz et al.l d2006l) and 
Mrk 817; iDennev et al.l (120101) ). In this case the mean spec¬ 
trum should be used instead. 

The quality and uncertainty in the width measurement 
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Figure 13. The predicted scatter in the recovered lags and recovered fraction for (from top to bottom) the long (2 extra years of 
observations), weekly (weekly cadence of spectroscopic measurements), full season (9 month observation season for the last 3 years of 
survey), year (hypothetical 5 year survey), goal (0.03 mag uncertainty in spectroscopic measurement), year+goal and weather (3-5 epoch 
of spectroscopic data were lost over the 5 year period) extensions compared to the baseline OzDES results (grey underlay, and top left). 
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Figure 14. The recovered fraction of the weekly [dotted; weekly 
cadence of spectroscopic measurements], full season [dashed; 9 
month observation season for the last 3 years of survey], long 
[dot dashed; 2 extra years of observations], and year [dot dot dot 
dashed; 5 year survey with no seasonal gaps and monthly spectro¬ 
scopic cadence] survey extensions as a function of the observed 
time lag compared to baseline survey [solid with grey shading; 
corresponding to the median value and inner 68 percentile values 
from bootstrap resampling]. 


are dependent on the l ine width characterisation used and 
the SNR of the spect ra (jPennev et al.l2009l ; iFine et ahll 20 id : 
I Jensen. Jens M . The common ways to characterise the 
velocity dispersion from the line width are the full width 
at half maximum (FWHM), the second moment of the line 
(crime; otherwise known as the li ne dispersio iEI ), and the 
inter-percentile velocity (IPV; IWhittla Il985h . We expect 
the fractional uncertainty in AV, measured from the mean 
spectrum to be <5% (< 0.04 dex in Mbh). However, the 
rms spectrum signal has an extra dependence on the level 
of variability, which will impact the S/N of the line pro¬ 
file and therefore the accuracy and precision of the line 
width measurement. Compared to the current reverbera¬ 
tion mapping sample, of primarily low luminosity objects 
(/t {log [ALa( 5100A)] } = 43.5), the level of variability in the 
OzDES sample is likely to be smaller. Therefore, we expect 
the line width uncertainties measured from the rms spectra 
of these quasars to be higher, with fractional uncertainties 
of order of 20-30% (0.15-0.20 dex in Mbh; K. D. Denney, 
private correspondence). Additionally, in some objects we 
expect the emission line flux errors per pixel to be on the 
same order as the variability signal, and consequently, the 
rms spectrum signal will be too low to resolve unambigu¬ 
ously. 

It should be noted, that care should be taken when mea¬ 
suring the mean spectrum to ensure that any non-variable 
components, like narrow emission lines, are modelled and re¬ 
moved. This is particularly important for Civ. It has been 


8 The ‘line dispersion’ is distinct from ‘velocity dispersion’ and 

is defined as rr 2 - /A\ 2 - /A 2 \ - S x2F W dx _ ( I *F(\)dX ' 2 
is defined as, cr line — (A) (A ) — j F( ^ dA [ / F(X)dX 
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Figure 15. The accuracy a a of, from top to bottom, the weekly 
(weekly cadence of spectroscopic measurements), full season (9 
month observation season for the last 3 years of survey), long (2 
extra years of observations), and year (5 year survey with no sea¬ 
sonal gaps and monthly spectroscopic cadence) survey extensions 
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Figure 16. Forecast of lag detections for the OzDES rever¬ 
beration mapping campaign. The different colours and symbols 
correspond to the different emission lines (black circles-H/3, red 
triangles-Mgn, blue crosses-C iv). 

found that using line dispersion is more robus t than FWHM 
in th i s case, as the pr ofile changes with width (jRafiee &; Halil 
HoH; d enne 3 Sm3). Additionally, when using the mean 
spectrum a different f-factor is required as the mean spec¬ 
trum width tends to differ from the rms spectrum width 
systematically. 

When absorption or contamination is present in the 
emission line profile, a rms spectrum can be created using 
individually modelled spectra (see Section [3.4.311 . However, 
if the SNR in each spectrum is not sufficient to accurately 
decompose the individual spectra, the mean spectrum is 
modelled instead and the width is measured a c cordin gly. 
This technique was performed bv lDe Rosa et all (2015) for 
NGC5548. 

On average, cfr/R ~20% (0.08 dex in Mbh) and no 
significant dependence on redshift is apparent. Thus, for 
the baseline OzDES setup we expect a median formal un¬ 
certainty in Mbh of 0.09 dex (0.16-0.21 dex) in a random 
sample of OzDES targets using the mean (rms) spectrum, 
ignoring the uncertainty in /. 

5.2 Recovery of the R — L relationship 

One of the major scientific goals of this survey is to derive 
the R—L relationship for all three emission lines. In this sec¬ 
tion, we investigate how to optimally select our target quasar 
sample to recover the most accurate and precise R—L rela¬ 
tionship. To recover the R—L relationship for all three lines 
we require a substantial calibration sample with lag mea¬ 
surements for two emission lines so that the relationships 
can be put on the same relative luminosity scale. This is 
crucial for both black hole mass es timates and any a ttempts 
to use quasars as standard candles dKing et al.l2014l ) . There¬ 
fore we decided to monitor all 69 of the quasars that fall into 
a redshift range where two lines can be observed simultane¬ 
ously (i.e. 0.54 < z < 0.62 and 1.78 < z < 1.96). 

The optimal strategy for constraining the R — L rela¬ 


tion is to observe objects over a broad luminosity range. 
To do so we need to observe both faint low -2 objects and 
bright high -2 objects. However, we also want to optimise the 
number of recovered lags. For the rest of this section, we de¬ 
termine whether selecting our targets randomly to cover a 
broad range of properties, or based on their expected accep¬ 
tance fraction or accuracy, leads to a better estimate of the 
R — L relationship. We separated our ten AGN mock cat¬ 
alogues into two groups, a training sample and a observed 
sample. The training sample was used to calculate which 
redshift-magnitude bins had the highest recovery rates and 
accuracy. AGN were then selected from the observed sam¬ 
ple according to their magnitude and redshift and used to 
calculate the R — L relationship for each emission line. The 
separation was made in an attempt to avoid any biases in 
the resulting constraints from particularly favourable or un¬ 
favourable bins that do not follow the general magnitude 
and redshift trends. 

The training sample was broken into separate redshift 
and magnitude bins and the recovered fraction and accuracy, 
(ja, was calculated for each bin. The resulting distributions 
were quite noisy, so to remove the influence of spurious bins 
that may skew our results we fit a low order polynomial 
surface to both the recovered fraction and cta training dis¬ 
tributions, using the ID 10 function SFIT. Due to the signif¬ 
icant differences between Mgn and Civ lags, we split the 
redshift range into a low- (z < 1.78) and a high-redshift 
(z > 1.96) group and fit each group separately. The recov¬ 
ered fraction distribution does not appear to have any major 
structure besides a decline towards bright high-redshift ob¬ 
jects, so we simply fit a linear distribution. The distribution 
of cta does exhibit several significant features due to the lag 
dependences found in the previous section, so we fit a 3rd 
order polynomial surface to the low redshift group and a 4th 
order polynomial surface to account for the apparent struc¬ 
ture in the high redshift group. The residual distributions of 
these fits did not show any significant underlying structure. 

The observed sample of 500 sources was then chosen 
from specific redshift and magnitude bins according to ei¬ 
ther recovered fraction or accuracy restrictions. The crite¬ 
ria tested were: (a) the recovered-percentage/acceptance was 
greater than 50%; (b) The recovered-percentage/acceptance 
was greater than 40%; (c) A was less that 0.05; or (d) A 
was less than 0.10. 

For each sample we fit a R — L relation of the form 

log(i?/l It year) = K + a\og{\L\/l erg s -1 ), (13) 

assuming the luminosity of each object is known exactly. The 
a and K parameters were then calculated for 1000 different 
possible samples of observed quasars. The resulting median 
estimates of a and K and thei r uncertainti e s are shown in 
Fig. 1171 and Table [3] Following iBentz et al.l (120131 ). we also 
define the scatter in the sample as the standard deviation 
of the residuals around the best fit relation, and its value is 
given in Table [3] Table [3] also includes the median observed 
and expected recovered number of H/3, Mgn, and Civ lags 
for each selection criterion. 

For the H/3 case, there appears to be a systematic bias 

9 Interactive Data Language (Exelis Visual Information Solu¬ 
tions, Boulder, Colorado) 
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towards a shallower slope in the recovered R — L relation¬ 
ship than expected. This is due to the poor lag accuracy 
for very short lags, associated with the lowest luminosity, 
low-redshift objects. We attempted to minimise this effect 
by restricting our sample to accepted grade 1 lags only. The 
resulting a and K values are shown by the triangles in Fig. 
El After this quality cut, the bias is no longer apparent and 
the recovered values are more tightly constrained. Therefore 
we should only include the objects with clear RM lag signals. 
We can also note, that the previous RM campaig ns constrain 
t he l ow luminosity, low redshift objects well dBentz et al.l 
1201311 and their contribution to the R — L relation fit has 
not been considered here. 

The Mg 11 R — L relation constraints for both the full ac¬ 
cepted sample and accepted grade 1 only sample were found 
to be reasonably consistent, although the scatter around the 
grade 1 best fit was ~30% smaller than for the full sam¬ 
ple. Therefore the Mgn sample is not systematically biased 
by spurious lag values, but the higher constraining power 
gained by including more objects is balanced by the increase 
in the scatter around the relationship. 

Like H/3, the recovered C IV R—L parameters also suf¬ 
fered from a bias towards a shallower R — L slope with the 
whole accepted sample, which is resolved by restricting the 
sample to grade 1. The resulting constraints accurately re¬ 
cover the input R — L relation parameters in all cases except 
for samples chosen for acceptance. This is due to the reduced 
number of CIV lags used for the parameter constraints. The 
original bias towards a shallower R — L slope is likely due 
to the short luminosity baseline for CIV R — L relationship, 
which is more sensitive to misidentified lags than the corre¬ 
sponding Mg 11 R — L relation fit. Despite both lines sharing 
a similar lag distribution and therefore overall lag accuracy. 

The accuracy and precision of the parameter estimates 
appear to be relatively independent of whether the sam¬ 
ple was selected randomly or with an accuracy constraint. 
However, a sample chosen for high recovered fraction created 
tighter constraints for H/3 and poor and possibly misleading 
constraints in Mg 11 and C IV. This is primarily driven by the 
number of lags recovered, although the luminosity baseline 
of the R—L relation fit is also restricted by this choice, which 
will affect the recovery of the R—L constraints. Therefore 
in order to accurately and precisely recover the R — L re¬ 
lationship for all three lines it is preferable to get an even 
distribution of targets over the total redshift and magnitude 
distribution rather than maximising the total number of re¬ 
covered lags. 


5.3 Extensions 

We also investigated the improvement in both ov/t and the 
recovered R—L parameters associated with the different 
survey extensions using a random target selection process 
(equivalent to the ‘None’ selection process in previous sec¬ 
tion). The results are summarised in Table [4] The loss of 
epochs due to weather had a universally detrimental effect 
on the recovery of Mbh and the R—L parameters. The pre¬ 
cision in the black hole mass measurements was most signifi¬ 
cantly improved by a reduction in the spectral measurement 
uncertainty (goal), however great improvement was also ob¬ 
served in the case where no seasonal gaps are present (year). 
The ‘goal’ extension is expected to have a median Mbh un¬ 


certainty of 0.05 dex, assuming no improvement in the AV 
measurement. Interestingly, the precision in mass from the 
‘goal’ case was found to be superior to the ‘year+goal’ case, 
though the difference is not highly significant. The differ¬ 
ence is likely due to the slightly longer photometric baseline 
of the baseline OzDES survey from the DES science verifica¬ 
tion monitoring in 2012. Note that in the previous section we 
found the accuracy of the lag measurement to be consistent 
between the ‘goal’ and ‘year+goal’ extensions. 

The R — L relationship parameter constraints for all 
lines were most improved from a reduction in the measure¬ 
ment uncertainty (goal; Fig.Q+J. The precision and accuracy 
in the H /3 and CIV R—L parameter constraints were also 
significantly improved in the cases where the seasonal gaps 
were reduced or removed completely (i.e. year and full sea¬ 
son). In the Mg 11 case, the distribution of constraints for the 
‘weekly’ extension exhibited a large tail towards a shallower 
R — L relationship. This is a consequence of the enhanced 
lag uncertainty created by the seasonal gap. If only grade 
1 lags are used in the construction of the R — L relation¬ 
ship, this tail disappears. Additionally, the R—L parameter 
constraints were only marginal improved by extending the 
survey by an extra two years (long). 


6 DISCUSSION AND CONCLUSIONS 


The baseline OzDES survey is expected to observe ~ 500 
AGN for the full five year period, spanning a redshift range 
of 0 < z < 4 and luminosity range of 10 39 < A L\ < 10 4 ' 
erg s _1 . If the final 500 AGN are randomly selected from 
the target quasar catalogue we expect a lag recovery rate 
of ~ 35 — 45% (Table 0. This would represent a four-fold 
increase in the number of measu red lags compared to the 
current sample dBentz et al.ll2013l 'l, and a more than ten-fold 
increase in redshift range. A higher acceptance rate (~ 60%) 
can be achieved if the final targets are selected based on their 
expected lag length and their current light curves. 

This sample will enable direct Mbh measurements over 
a broad range of quasar properties, environments and black 
hole masses. It will also help constrain the R — L relation¬ 
ship for multiple emission lines and test the robustness of 
this relationship over a broad AGN population, including 


an investigation into the recently observed Eddin gton ratio 
depe ndence of the R—L relationship discovered by (iDu et al.l 
120151) . The current reverberation mapping sample is biased 
towards local, low luminosit y objects, which are not repre¬ 
sentative of typical quasars ( Richards et al.l l2011h This had 
lead to concerns that the existing RM sample may be biased 
compared to the broader quasar populati on (e.g., Shen et ahl 
120081 ; [Richards et al.ll201ll : D enne 12012IL The OzDES sam¬ 
ple will probe a similar redshift and luminosity range to the 
quasar samples in which the secondary mass estimate tech¬ 
niques are applied and thus minimise any potential biases. 

The baseline OzDES reverberation mapping campaign 
is predicted to accurately recover the R — L relationship 
for all three lines when only the accepted grade 1 sample is 
used and a wide distribution of targets are selected (Figure 
El)- In general, the spread in the recovered R — L param¬ 
eters is larger than the uncertainties associated with the 
input relationships. However, the current Mgn R—L rela¬ 
tionship is not based on direct Mgn lag measurements (see 
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Figure 17. The median recovered gradient (a: left) and intercept right) of the R — L relation for H/3 (top), Mgn (middle) and Civ 
(bottom) for both the full accepted sample (circles) and the accepted - grade 1 sample (triangles) for each sample selection method. The 
horizontal line and grey shaded regions show the input values for a and K and their current uncertainties. 










































































Mgll 


20 A. L. King et al. 




Figure 18. The median recovered gradients (u: left) and intercepts right) of the R — L relation for the various survey extensions. 
Results are shown for grade-1 H/3 (top), Mgll (middle) and Civ (bottom) lines. The horizontal line and grey shaded regions show the 
input values for a and K and their current uncertainty. 
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Table 3. Recovered R — L parameters for the five selection criteria. 


Selection Criteria 

Line 

K 

a 

Scatter 

Total lags 

Accepted lags 

None 



H/3 


U.4o_o 04 

pj 1 91 +U.U^U 
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80±8 

51 ± 6 
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H/3-grade 1 
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0 52"*" 0 ' 04 

u.O^_o.03 

0 107+ 0 036 
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> 

40% 

H/3 

-18.31?'? 

o.^O-o 01 

u.iz/_ o 010 

266±10 

169 ± 10 
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> 
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95 ± 8 
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> 
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0.1081!;°!® 

91±5 

37 ±6 

Acceptance 

> 

50% 

Mgll-grade 1 

-26.41?;® 

0 f}Q+ 0 - 08 
o.DO-ono 

0.0811!;°!! 

91±5 

17 ±5 

Acceptance 

> 

50% 

Civ 

-22.0lg® 

0 52”*” 0 ' 13 

O.OZ-o 15 

0 1 55+ 0 - 027 
0.10O—0 027 

62±0 

26 ± 3 

Acceptance 

> 

50% 

C iv-grade 1 

oq 1 +8.0 

OO.l-n 9 

pi '77’+0.26 
u -' '-0.18 

0 1Q7+0 030 
0-10' -0.042 

62±0 

11 ±2 

A < 0.05 



H/3 

-i8.4ii;i 

0 4 [ 5+ 0 - 05 
O.4O-0 04 

0 11R+0- 016 

U.ll D _o.oi7 

75±8 

49 ± 6 

A < 0.05 



H/3 - grade 1 

-21.3li g 

0 52 +0 - 04 
O.o^_o.o4 

pj 1 1 q+ 0.022 

u.iiO-o.022 

75±8 

27 ±5 

A < 0.05 



Mg ii 

-26.0ll;l 

0 62"*" 0 08 
U.OZ-o 04 

o.nil!;!!! 

241±11 

100 ± 10 

A < 0.05 



Mgll-grade 1 

-26.31?;! 

0 «q + 0 - 04 
U.OO-0.05 

o.oselg;!!! 

241±11 

41 ±6 

A < 0.05 



Civ 

-20.611;, 

0 4q+ 9 - 93 

0-40-0.05 

n 199+0.020 
U.1ZZ-0 016 

253±10 

78 ± 9 

A < 0.05 



C iv-grade 1 

99 9+2.7 

0 53+ 0 06 
O-OO—0.06 

0.098lg;°!5 

253±10 

28 ± 5 

A < 0.10 



H/3 

-i8.8!£5 

0 46"*” 0 04 

0.40-0 04 

p) 114+0- 024 

0.114_o 018 

69±7 

44 ±6 

A < 0.10 



H/3 - grade 1 

—21.311! 

0 52 - * -0 ' 04 
O.o^_o.o3 

0 ioi+o - 040 

U.1U1-0 022 

69±7 

25 ± 5 

A < 0.10 



Mg ii 

-25.4120 

0 61 + 0 - 10 

0.01_o.05 

0.1191!;°!° 

212±11 

87 ±8 

A < 0.10 



Mgll-grade 1 

-25.8ll;l 

0 62"*" 0 ' 06 
u - dz - 0.07 

0.092lg;°!° 

212±11 

35 ± 5 

A < 0.10 



Civ 

-21.2II5 

0 5Q+ 0 - 04 
U.OU-0.04 

0.1301!;°!® 

287±11 

98 ±9 

A < 0.10 



C iv-grade 1 

-22.61®;° 

0 54+ 0 - 06 

O.O4-0 06 

0.0981!;°!® 

287±11 

34 ±5 


Table 4. Recovered R — L parameters for the survey extensions. 


Extension Rb .0 ckh/? Scattern/3 A/Mgii ^Mgii Scattering n Rciv a:civ Scatterc iv ^r/R 


Default 

-19.61!;! 

10 

0 9 

6 d 
+ 1 

00 

0 

n 199+0.020 

o.izz-o Q17 

-25.6111 

0 ni+0- 08 
u - D1 -0.05 

pj 190+O.OI8 

u.±zo_0 016 

-21.il!;® 

P) cpj+0.04 
U.OU-004 

0 I 29”*~ 0013 
U.1ZU-0 014 

P) 91 Q + 0127 
U.Z1O-0 10 4 

Weather 

99 o+4.3 
^■°-6.5 

0 54+ 0,15 

W.04_o_ 10 

pi 1 ro+0.050 
o.J-OO-o 039 

-25.1I?;? 

0 60" 1 ” 0 04 
u.ou-o 05 

0.1081!;!!! 

-20.61!;! 

0 49 - *" 0,03 
u.^-0.04 

0 107+ 0 - 009 
U.1U/ _o 009 

0 211" 1 ” 0 ' 124 
u.zii-o in 

Goal 

-2l.5lg;| 

f) co+0.02 

U.OO-o 02 

0 Oh'S - * -0 023 
U.UOO_o 025 

—25.5l!J 

0 61 +0 01 
u - ul _0.oi 

0.046lg;!!! 

-23.3lg;J 

P) cc+0.01 
u - oo _0.01 

0 046 +0 005 
u.u^u_o 005 

0 07 c 5+°' 063 

u.u/o_o 030 

Long 

-19-711:3 

0 48”*” 0 ' 05 
^•40-0.05 

0.167lg;g£ 

-23.71!;! 

0.57lg;!! 

0 IO8”*” 0 015 
u.iuo_o 013 

-21.81!;? 

0 52 - *" 0 02 
u - oz -0.03 

0 097 - * -0,008 

u.uy < _o 008 

0 ini5+ 0 - 084 

u.ioo_o 073 

Weekly 

-21.61!° 

0 53+ 0 03 
u.oO-o.02 

0.0921!;!!° 

-24.51°;! 

0.591!;°? 

0 091 ”*” 0 027 
u.uui-o 016 

-20.81!;! 

P) cn+0.03 

U.OU_oo4 

o.ioeig;!!! 

0 119+ 0089 
u.iiy_ 0 .045 

Full Season 

-21.9I?;! 

0 53+ 0 - 03 

u.oo_o 02 

pj 1 9^+0.032 
U.izO-o 027 

-25.21!;! 

0 60"*” 0 01 
u ‘ ou -0.01 

0 1 1 c+0.044 
U.11O-0 049 

—24 l - * -0,8 

^4.i_o 8 

0 57+ 0 02 
u -° * —0.02 

0.0821!°!? 

pj 190+O.O86 
u.J-^O-o 055 

Year 

-20.6lg° 

0 5Q+ 0 - 01 

u.DU-o 01 

0.0581!;!!? 

—24 Q +1 - 5 

^4.U_o 9 

0 60”*” 0 02 
u.ou_o 03 

0 04Q+ 0 007 

u.u^y_o 005 

99 7+O.7 

ZZ, '-0.7 

0 54+ 0 - 02 
u -°4_o .02 

0 060 +0 ' 008 
u.uou_o 008 

0 Q77+ 0 - 057 
u ' u ' ‘ -0.035 

Year+goal 

-20.91!;? 

Q tri+0.02 

u - 0i -o.oi 

0.0731!;!!® 

-25.91!;! 

p) fi2”*” 0,02 

u.oz-o 02 

0 nnQ+0 046 
U.UOo —0.013 

-22.91!;? 

0 ^4+ 0 - 01 
U * °4-0.01 

0 046 ”*” 0 ' 006 
U.U4D_o 005 

0 Q7«+ 0 060 
U.U *0—0.034 


Note: In this case, the H(3 sample corresponds to the H/3 - grade 1 sample. 
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iTrakhtenbrot fe Netzeijl2012lh while OzDES will constrain 
it directly, and the C IV R — L relationship is based only on 
a small number of objects. 

The large spread in the C IV R — L relationship parame¬ 
ters is due to the combination of the low acceptance rate, low 
lag accuracy and relatively short luminosity baseline of the 
C IV sample. However, any of the survey extensions signifi¬ 
cantly improve the R — L relationship parameter constraints 
for all lines (Figure [T5] and Table [3J. Accurate estimates of 
all three R—L relationships are crucial for single epoch mass 
estimates and measuring distances to the highest redshifts. 
Consequently, there are significant gains from pursuing one 
or more of the survey extensions we have simulated. 

One of the greatest concerns with using AGN as stan¬ 
dard candles is whether the R—L relationship evolves in 
redshift. The OzDES sample may enable an investigation 
into any trends in the R — L due to redshift, metallicity, 
Eddington ratio and many other properties. If the R — L re¬ 
lationship appears to be consistent over the observed sam¬ 
ple of quasars, then the OzDES reverberation sample will 
provide the first physically motivated distance measurement 
based on a single method from the present day back to red¬ 
shift four. Unfortunately, the statistical power of the OzDES 
sample is not expected to be competitive with existing cos¬ 
mological probes as the predicted uncertainty of the R — L 
will still be too large to rival the precision in current SNe 
and BAO measurements. However, it has the potential to 
uncover unexpected expansion behaviour if large deviations 
to ACDM are present at high redshifts and will provide a 
strong base for future surveys. 

The Sloan Digital Sky Su rvey (SP SS) is cu rrently run¬ 
ning a comparable campaign (IShen et al.ll2015f l on a much 
shorter time scale (6 months in 2014). They are observing 
849 quasars in a 7 deg 2 field of view using the SDSS-III 
Baryon Oscillation Spectroscopic Survey (BOSS) spectro¬ 
graph. Their sample is flux-limited to i ps f = 21.7 mag (best- 
fit point-spread function (PSF) magnitudes from SDSS), has 
~30 epochs of spectroscopic data over the duration of the 
survey (an ~4 day cadence), photometric monitoring ap¬ 
proximately every 2 days, and includes quasars up to red¬ 
shift 4.5. They expect to recover lags for 10% of their sample 
out to a redshift of 2, with a possible extension to z~4 with 
the inclusion of 3 years of photometry obtained over 2011- 
2013 from the Panoramic Survey Teles cope an d Rap i d Re - 
sponse System 1 Survey fPanSTARRS. iKaiser et al]|2010l) . 
The OzDES and SDSS reverberation mapping campaigns 
will produce quite complementary measurements. SDSS has 
finer temporal sampling and a shorter timeline, enabling 
the recovery of a broad range of faint AGN with shorter 
lags, while the DES/OzDES sample will be able to more ef¬ 
ficiently recover the brighter and higher redshift AGN with 
longer lags. 

6.1 Target selection criteria 

Selecting our sample based on the expected accuracy did 
not significantly change the predicted R—L parameter con¬ 
straints, although, when the sample was chosen for a higher 
recovered fraction we observe a significant tightening in the 
H p constraints and degradation of the Mgn and Civ con¬ 
straints (Figure HT1) . This is primarily due to the relative 
number of lags used to constrain the R — L relationships. 


We find it very advantageous to prioritise targets with mul¬ 
tiple lines present in their spectra to calibrate between the 
H j3 R—L relationship, which is quite tight ly constrained by 
the current RM sample dBentz et al.ll2013l) . with the R — L 
relationship of the higher ionisation lines. It is also impor¬ 
tant to cover a diverse range in redshift and magnitude. 

6.2 Extensions 

Of all the survey extensions, improving the measurement 
uncertainty, through better spectroscopic calibrations, was 
the most efficient means of improving the overall results. 
This is because the constraints on each emission line data 
point are much stronger, leading to less ambiguity in the 
lags. The goal measurement uncertainty of 3% is quite op¬ 
timistic but may be possible to achieve. During the span of 
the survey so far, major upgrades have been made in both 
the AAOmega instrument 0 and its pipeline Another 
method of reducing these uncertainties is to increase the 
number of calibrating F stars monitored in each field. 

In terms of improved scientific results, reducing the 
measurement uncertainty was closely followed by closing the 
seasonal gaps, even by 3 months. By extending the observing 
season beyond the 6-month season of the standard survey we 
break the ‘half year degeneracy’ and allow common variabil¬ 
ity features to be probed by both continuum and emission 
line light curves. However, the improvement in the lag re¬ 
covery was generally limited to the sh orter lags. This re - 
sult is consistent with previous findings dHorne et alj|2004f l. 
namely that the recovery of longer lags requires a longer 
baseline of observation rather than finer sampling. 

Extending the observation timeline by 2 years (long), 
only marginally improves the predicted scientific results, de¬ 
spite expectations to the contrary. A longer program does 
enable longer lags to be recovered (as seen in Figure El), 
which allows a broader luminosity baseline from which to 
constrain both the Mgn and Civ R—L relationships. How¬ 
ever, due to the random selection process we employed and 
the lower number density of brighter objects, this extension 
had a trivial influence on the R — L parameter constraints. 
It is likely that brighter objects will have relatively higher 
priority in our target selection (Fig. [T]), so in the real survey 
this improvement may be more substantial. The longest lag, 
highest redshift objects are also best studied with a long 
program like OzDES, and thus there are clear advantages 
to monitoring as many of those objects as possible in the 
hope that an extension of the survey will prove possible. 
The overall ranking of the extensions is shown in Table [5] 
The rank is determined from the median value of each cri¬ 
terion and the final ranking is calculated from the overall 
sum of the other ranking values. The criteria tested were, 
the number of additional hours of observation required (re¬ 
ferred to as resources), recovered fraction, accuracy (<ta), 
precision in lag measurement (oy/r), and the precision and 
accuracy of R—L parameter recovery. The extensions from 
best to worst ranking are: Goal, Year, Full Season, Long and 
Weekly. The final ranking was based on the overall sum of 
the other ranking values. 

10 http://www.aao.gov.au/science/instruments/current/status 

11 http://www.aao.gov.au/science/software/2dfdr 












A GN Reverberation Mapping Project 23 


Table 5. Survey extension rankings. 


Extension 

Resources 
(Add. Hours) 

Recovered 

Fraction 

<?A 

ov/t 

R — L Accuracy 

R — L Residual Scatter 

Overall Rank 

Goal 

1 

1 

1 

1 

1 

1 

1 

Full season 

2 

3 

2 

4 

3 

3 

3 

Long 

3 

4 

5 

5 

4 

4 

4 

Weekly 

5 

5 

4 

3 

5 

5 

5 

Year 

4 

2 

3 

2 

2 

2 

2 


The rankings (1-best to 5-worst) are given for the required number of additional hours of observation (resources), recovered fraction, 
accuracy, ov/t, and R — L parameter recovery. The rank is determined from the median value of each criterion and the final ranking is 

calculated from the overall sum of the other ranking values. 


Losing an additional 3-5 scheduled epochs of spectro¬ 
scopic data did not significantly affect the recovered fraction 
and accuracy of the overall mock catalogue sample, but did 
severely affect both the precision of black hole mass mea¬ 
surements and R — L relationship parameters. Therefore it 
is important to minimise the number of epochs lost over the 
observation period. One way OzDES is working to minimise 
potential losses is by working in close collaboration with the 
2dFLenS survey, also using the 2dF instrument, to make 
supplementary observations of the SNe fields when weather 
restricts OzDES observations. In exchange, OzDES will ob¬ 
serve 2dFLenS fields when the DES SNe fields are at high 
airmass. 


6.3 Alternative analysis techniques 


On completion, OzDES will be one of the longest running 
reverberation mapping campaign to date and in number of 
AGN monitored, it is second only to SDSS. Nevertheless, the 
expected number of OzDES spectroscopic epochs is small 
com pared to traditional r everberation mapp i ng campaigns 
(e.g. j Peterson et al.||2002l: iBentz et al.ll2009al ; iDennev et al.l 
120091 ; Barth et al. 2011 ). This leads to the relatively low re¬ 
covered fraction and accuracy predicted. However, one way 
to maximise the output of the OzDES data is to stack the lag 
signals for multiple objects of similar reds hifts and magni¬ 
tudes (‘composite rever b eration mapping’, iFine et ahl 120121 : 
iBrewer fe Elliottl I20l4l . IFine et ahl (120121 ) found that by 
stacking the continuum and emission-line light curve cross 
correlation signals of objects with similar redshifts and mag¬ 
nitudes, and therefore similar expected lags, a mean signal 
is recovered even if no lag signal is present in the individ¬ 
ual cross correlations. The large number of OzDES targets 
makes it a good sample to perform this type of analysis. 

We will also perform reverberation mapping on other 
lines in the spectrum beyond H/3 A4861, Mgn A2798, and 
C IV A1549. Only these three lines were mentioned in this 
work as they are the three lines traditionally used for single 
epoch mass determination, and calibration of the R-L rela¬ 
tionship for single epoch masses is one of the main science 
drivers of this survey. 


6.4 Limitations of survey simulations 

There are several limitations to our survey simulation setup. 
The first is our continuum luminosity determination for the 
individual objects. We have simply used the existing SDSS 
template, and have not taken into account host galaxy con¬ 
tamination or extinction. Our choice of quasar template was 


based on the similar redshift and magnitude range covered 
by the SDSS sample and the OzDES target sample. The 
slope of the spectral energy distribution (SED) has been 
found to vary consid erably between indiv idual objects and 
different samples (e.g lRichards et al.ll2006l ) so using a single 
template is a simplification. However, it is sufficient for our 
use as we are investigating the efficiency of the survey for 
the bulk of the quasar sample. 

Not taking into account host galaxy contamination 
overestimates the bolometric luminosity of the AGN, result¬ 
ing in reduced sensitivity in our observations, an overesti¬ 
mation of the lag length, and an underestimation of charac¬ 
teristic variability. However, this is only significant for AGN 
of similar luminosities to their host, and in general, the ex¬ 
pected decrease in lag recovery due to lower sensitivity in 
light curve variation measurements is counteracted by the 
increased variabili ty and sho r tened lag length of lower lu¬ 
minosity quasars. IShen et al.l (120151 ) took this into account 
in the SDSS reverberation mapping campaign by assigning 
a constant host galaxy contribution of 8 x 10 43 ergs -1 at 
5100A and negligible contrib ution at L %nnn i aid L 1350 ^. 
If we follow the approach of [Richards et al.l J2006I) (with 
LbotJR Edd = 0.25 based on the results of iKollmeier et al.l 
iSi)) to estimate the host galaxy contribution of our sam¬ 
ple, we found that only low redshift objects are expected to 
have a significant host galaxy component and the majority 
of our targets will not be affected appreciably by neglecting 
host galaxy light due to their intrinsic luminosities. To con¬ 
firm this, we tested how the recovered fraction and accuracy 
of a mock sample were affected by host galaxy contamina¬ 
tion and found the results to be consistent with the baseline 
simulation. However, the inclusion of host galaxy contami¬ 
nation does result in a systematic drop in the expected lag 
length and again the effect is dominant at low-redshifts. 

Ignoring internal extinction also underestimates the 
bolometric luminosity. However, this is only expected to be 
a 20% luminosity correction base d on extinction es timates 
for the SDSS DR9 quasar sample (I Paris et al ]Eml) , which 
will not cause significant changes to the simulation results. 

Another potential limitation of our simulation is our 
choice of transfer function. A top hat transfer function is a 
good conservative choice in this type of investigation, as it 
spreads the lag response more dramatically than more Gaus¬ 
sian transfer functions. Nonetheless, it could be argued that 
by choosing a top hat transfer function, when JAVELIN is 
based on a top hat transfer function, we are biasing our 
results. Our choice of top hat width may also affect our 
results. To test both of these issues we consider two alter¬ 
native transfer functions. First, we considered a Gaussian 
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transfer function, with a width of O.l r, motivated by ex¬ 
isting reverberation mapping data ll Gi ier et~al IB) and 
the SDSS survey simulations (jShen et al.ll201fll h The recov¬ 
ered fraction and accuracy were consistent or slightly better 
than the top hat predictions fFig. 1191) . We also investigated 
a worst-case top hat transfer function scenario, where the 
top hat width is 2r instead of 0.2r. In this case we do see 
a significant drop in the recovered fraction and accuracy 
due to the extreme smoothing of the emission line signal. 
This case can be considered as the worst-case result for the 
baseline OzDES survey. We expect negligible effects from 
changing the width of the t ransfer func tion by a factor of 
2-3, based on the findings of IS hen et al.l d2015l) . who found 
that lag recovery is not significantly affected by changes in 
this range. 

Additionally, the choice of top-hat scaling factor will 
have a strong effect on the lag recovery. Our ch oice of unity 
was based on lZu. Kochanek fc Petersonl (1201 il l findings for 
NGC5548, which were on the order of one. However, they 
also found a possible correlation between luminosity and am¬ 
plitude that we have neglected in our analysis. The choice 
of amplitude directly relates to the responsivity of the line. 
As mentioned previously, we theoretically expect different 
emission lines to respond b y different degrees to variations 
in the ionising continuum ( [Goad. O’Brien fe Gondhalekail 

1 19931 ; iKorista fe Goadll2000l . l2004lf and this behaviour is wit- 

nessed observational ly with sometimes contradictory trends. 
ICackett et al.l d2015l l recently found that the responsivity of 
Mg II was low in NGC5548, and displayed virtually no re¬ 
sponse to the continuum variability over an period of 170 
days, despite a near-UV continuum variability amplitude of 
Fvar = 0.3 3 . This finding agrees with the observations of 
iGlavel et al.l ( 1991 1 who found a smaller variability ampli¬ 
tude in Mg II flux c ompared to t he other UV lines observed. 
On the other-hand. IWod d2008l l found low Mgn responsiv¬ 
ity in only one of their objects and reasonable Mg II respon¬ 
sivity in the other four objects. Therefore, defining the re¬ 
sponsivity for any line or object appears to be complex. 
Additionally, some objects have also displayed very non- 
linear responses to continuum variations (e.g., NG C7469 
IPeterson et, al.l [20141 : NGC5548 iDe Rosa et alll2015l durin; 
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the s econd half of the campaign; and J080131 _ 

1201511 . Until the mechanisms that drive changes in respon¬ 
sivity are well understood we can only simulate the transfer 
function in a reasonable fashion, as we have done in this 
work. Based on our knowledge from previous RM programs 
and the diversity of the OzDES survey, it is hard to deter¬ 
mine whether the OzDES will perform better or worse than 
predicted in these simulations. 

Although host galaxy contamination is not expected to 
significantly affect the recovered fraction or lag accuracy, 
not accounting for host galaxy contamination will affect 
the determination of the R — L relation parameters. Host 
galaxy contamination can be estimated, for the observed 
OzD ES data, using a combin ation of spectral decomposi¬ 
tion (|Vanden Berk et al.ll200fi h. and high r esolution imaging 


as performed by 


Bentz et al.l ( 2009al . l2013l l. 


6.5 Other considerations 

We expect certain sections of the extracted spectra to be 
suboptimal for flux measurements, in particular the region 
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Figure 19. The effects of using a very broad top hat (w = r) or 
a Gaussian (a = O.lr) transfer function on the simulated results. 
The results for the default O.lr top hat transfer function distri¬ 
bution is shown in grey. For comparison, the baseline scenario 
statistics are: Accepted: 39.7 ± 0.7%, Misidentified : 6.4 ± 0.6%, 
and <ta = 0.083 ± 0.004. 


near the dichroic split of the AAOmega spectrograph at 
570nm and the Fraunhofer A+B absorption bands at 759nm 
and 687nm. We recommend prioritising targets whose emis¬ 
sion lines avoid these wavelengths by 2x the FWHM of the 
line. This corresponds to approximate redshifts of 2 ~ 0.17 
and a ~ 0.56 for H/3, 0 ~ 1.0 and 2 ~ 1.7 for Mgn and 
2 ~ 2.7 and 2 ~ 3.9 for Civ. Avoiding the Fraunhofer bands 
is of less concern than the dichroic split. 


6.6 Summary 

We have generated mock catalogues of AGN and correspond¬ 
ing light curves according to the expected OzDES sampling 
and quasar properties. We attempted to recover the input 
lag from the simulated light curves to quantify the efficiency 
and accuracy of the lag recovery. These results were used to 
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predict the expected performance and scientific output of 
the OzDES reverberation mapping project and several pro¬ 
posed survey extensions. We expect OzDES to yield lags for 
~ 35 — 45% of the monitored quasars. The resulting direct 
Mbh measurements are expected to have formal uncertain¬ 
ties of 0.16-0.21 dex and the baseline OzDES reverberation 
mapping campaign will accurately recover the R — L rela¬ 
tionship parameters for H/3, Mgn, and Civ. However, sub¬ 
stantial improvements can be gained if we either increase the 
spectroscopic cadence, extend the survey season, or improve 
the spectroscopic measurement accuracy of the survey. 
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APPENDIX A: QUALITY OF OZDES SPECTRA 

To quantify the quality of the OzDES spectra we estimated 
the SNR of the emission-line flux in the rest-frame spectrum, 
as follows: 

(i) For each line, the continuum flux density ( S c ) is de¬ 
fined by a linear function of wavelength anchored by the 
median pixel value in the upper and lower continuum re¬ 
gions, defined in Table IAT1 

(ii) The line flux, Fl, is measured as the integrated flux 
above the linear continuum flux within the specified line 
integration wavelength limits such that 

F l = Y^[Si- S?], (Al) 

»i 

where Si is the uncalibrated flux density value at pixel i, and 
i\ and 12 are the pixel values that correspond to the defined 
minimum and maximum emission-line wavelength range. 

(iii) The associated variance in the line flux, g 2 {Fl), is 
measured as the integrated variance within the emission-line 
wavelength region. The SNR is then calculated as, 

SNR = —(A2) 

An example of a spectrum and the associated line sum¬ 
mation limits and continuum fits is given in Fig. lATl The 
SNR. for a single object can vary markedly between lines 
and epochs due to changes in observation conditions and in¬ 
dividual line strength. Therefore, we allocate an SNR value 
for each QSO based on the median SNR of the best mea¬ 
sured line in the spectra. Fig. lA2l shows the magnitude and 
redshift distribution of quasars with different SNR value cut¬ 
offs. There appears to be no marked difference between the 


Table Al. Line Summation Limits 


Emission 

line 

Line 

integration 
limits (A) 

Continuum 
window 
lower (A) 

Continuum 
window 
upper (A) 

H/3 A4861 

4810-4940 

4770-4800 

5100-5130 

Mg ii A2798_^ 

2700-2900 

2660-2700 

2920-2960 

C iv A1549 

1470-1620 

1440-1470 

1700-1730 



2600 2700 2800 2900 3000 

Wavelength (A) 


Figure Al. An example of the line and continuum decompo¬ 
sition implemented in the SNR measurement process. The dia¬ 
grams shows only the Mg II A2798 emission line. The spectrum is 
shown by the black solid curve. The associated flux density uncer¬ 
tainty, -y/Vais shown by the surrounding red shaded region. 
The line summation limits and associated continuum regions are 
shown by the dashed and dotted vertical lines, respectively. The 
linear continuum fit is shown by the solid line and the resulting 
SNR in the emission-line flux measurement for this example is 47. 


redshift distributions of high and low SNR objects; however, 
a significant shift in magnitude is observed such that higher 
SNR objects tend to be brighter, as expected. 

These signal-to-noise measurements are preliminary, 
and only serve to quickly quantify the expected distribu¬ 
tion of emission line measurement quality within the OzDES 
sample. Line strength measurements will improve with the 
more thorough methods we will utilize for the final RM anal¬ 
ysis, including decomposition of the spectra to remove con¬ 
tamination. 


12 The regions around Mgn are severely affected by Fe II con¬ 
tamination so this choice of pseudo-continuum wavele ngth range 
should be used critically (IVestergaard &; WilkesIbOPll) . 
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Figure A2. The top panel shows the distribution in characteristic emission-line flux measurement signal-to-noise (SNR) values for each 
object, the dots represent SNR>10 and crosses represent SNR<10. To achieve the baseline emission-line uncertainty of 0.1 mag a bare 
minimum of SNR>10 is required. The histograms show the magnitude distribution (left) and redshift distribution (right) of the whole 
sample (solid) compared to objects with SNR values greater that 10 (dotted), 20 (dashed), or 30 (dot-dashed). 




























